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Abstract 

We examine a wide class of multi-field inflationary models based on fields that decay or stabilize during 
inflation in a staggered fashion. The fields driving assisted inflation arc on flat, short stretches, before they 
encounter a sharp drop; whenever a field encounters such a drop due to its slow roll evolution, its energy 
is transferred to other degrees of freedom, i.e. radiation. The rate at which fields decay is determined 
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dynamically and it is not a free parameter in this class of models. To compute observables, we generalize the 



\ analytic framework of staggered inflation, allowing for more general initial conditions and varying potentials. 

By searching for generic situations arising on the landscape, we arrive at a setup involving linear or hilltop 
potentials and evenly spread out initial field values. This scenario is not more fine tuned than large-field 
models, despite the fact that many more degrees of freedom are involved. Further, the ^-problem can be 
alleviated. 

The additional decrease of the potential energy caused by the decay of fields provides leading order 

o ■ 

contribution to observables, such as the scalar and tensor spectral index or the tensor to scalar ratio, for 



which we derive general expressions. We compare the predictions with WMAP5 constraints and find that 



hilltop potentials are borderline ruled out at the 2cr-level, while linear potentials are in excellent agreement 
with observations. We further comment on additional sources of gravitational waves and non-Gaussianitics 
that could serve as a smoking gun for staggered inflation. 
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I. INTRODUCTION 



The embedding of single-field inflation within string theory is an active field of study, see 



for recent reviews. One of the most successful setups yet is the KKLMMT proposal 
5J]. A complete implementation is not easily achieved since many conceptual challenges need to be 
overcome. Among these challenges is the requirement that all moduli fields need to be stabilized, 
that the inflaton potential ought to be extremely flat and that corrections to the potential must be 
understood well enough in order to guarantee that the flatness is not spoiled once a field traverses a 
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super-Planckian distance. Remarakably, it seems possible to addresses these challenges within the 



KKLMMT proposal, but at the cost of fine-tuning and a construction requiring delicately arranged 
ingredients [g, 7| (see however [8|]). 

But if we turn our attention to a generic location on the landscape, single-field inflation appears 
unlikely: moduli-fields can be dynamical, potentials are generically too steep and they change 
significantly if the fields roll far enough. A way to circumvent some of these problems is to 



implement assisted inflation 



lfl 
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1 12 . [l3j, a type of multi-field inflation, whereby many 
fields are dynamical and assist each other in driving the inflationary phase. In this scenario, the 
potentials can be steeper without spoiling slow roll, because a given field sees the driving force of 
its own potential, but it is slowed down by the combined Hubble friction of all fields. In addition, 
if the number of fields is large, none of them needs to traverse a super-Planckian stretch in fields 
space, thus alleviating the ry-problem to some extent [9, [ill . [l2^ . From this point of view, multi- 
field models are more desirable and natural than single-field ones. Some recent realizations of 
assisted inflation are A^-flation Q, (fields are identified with axions arising from some KKLT 



compactification [16j of type IIB string theory, see also 



13, 



18fl), M5-brane inflation (inflatons 



are identified with distances between branes that are spread out over an orbifold) or inflation from 



tachyons 
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Of course, all moduli-fields need to be stabilized once inflation is over, that is moduli stabilization 
has to occur during or towards the end of assisted inflation. Current models of multi-field inflation 
are often described by an effective single-field model, see i.e. 22] for a review; it is further assumed 



that all fields stabilize coherently, leading to a sudden end of inflation and (p)reheating just as in 
the single- field case . However, the use of an effective single field is problematic for several reasons; 
first of all, there is no a priori reason to expect that fields decay (stabilize) all at once; indeed, 
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25, 



261 ] . where a coherent 



there are cases, such as in inflation from multiple M5-Branes 
decay (stabilization) is impossible 2 . A proper computational framework has to incorporate this 



feature, such as in 



27 



2a | where analytic tools for staggered inflation were developed (see also 



the numerical treatment in 



251 ]. where this effect was named cascade inflation). Another problem 



of an effective single-field description is its inapplicability during reheating: even if inflation ends 
coherently for all fields^they generically de-phase during reheating, resulting in the suppression of 



many resonances 



29, 



30|, albeit tachyonic preheating 



3J] might still work [30|, |3J, |3J] . Further, 



1 Note that reheating predominantly hidden sectors can be a problem for multi-filed models [23q . 

2 During inflation the M5-branes separate from each other until they encounter a boundary brane into which they 
dissolve, one after the other. 
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the production of isocurvature perturbations and non-Gaussianities cannot be computed correctly. 



To incorporate that fields decay or stabilize during inflation, we extend the formalism of [271 ] . 
employing again the coarse grained description proposed therein: by smoothing out the number 
of fields M(t) we can introduce a continuous decay rate T(t) = —J\f/N. This rate causes energy 
transfer from the inflaton sector to an additional component of the energy momentum tensor Tr , 
i.e. radiation. To avoid spoiling inflation, the ratio e = 3(f + w r )p r / (2p to t a i) ~ T/(2H) needs 



to be small. The resulting setup shares some similarities to warm inflation 
inflation 



34, 



351 ] or dissipative 



361 ] . without many of its shortcomings (see i.e. [37f] for some problems within warm 
inflation related to the additional frition term, which we avoid). A study of scalar perturbations 
within staggered inflation reveals that e appears alongside the usual slow roll parameters e and r\ 
in observables Ql. 
In 



271 ] several simplifying assumptions were made, such as identical initial field values and 
potentials for all inflatons. In this paper, we relax these conditions to allow for more general cases, 
but we still assume flat potentials. To be precise, we assume that the change of the potential 
energy in a given field over its slow roll phase is small compared to its potential energy. Within 
this setup, we compute the spectrum of gravitational waves Vt, the spectral index tit and the 
tensor to scalar ratio r. 

Equipped with this improved formalism, we construct new models of multi-field inflation moti- 
vated by generic potentials on the landscape [38], with J\f ~ I0 3 fields and e 3> e 3 : we arrive at 
models where fields can participate in assisted inflation until they encounter a sharp drop, at which 
point they decay or stabilize quickly. Thus, the decay rate is determined by the dynamics and it is 
not a free parameter. We assume similar potentials (hilltop or linear) for all fields, but allow for a 
narrow spread of masses or slopes. These types of potentials are motivated by expanding around 
either flat stretches or maxima on the landscape. Regarding initial conditions, we spread the fields 
evenly over the allowed interval and argue that this choice is the least fine tuned one. 

The scenario is similar in spirit to chain inflation 



39, 



40 



41 



42 



43, 



44 



451 ] . which is supposed 



to operate via successive, rapid tunnelling between different vacua on the landscape. However, 
instead of a rapid succession of first order phase transitions, which is difficult to achieve and can 
cause fatal problems in chain inflation related to bubble nucleation [44]], fields follow a smooth 
path through the landscape. The presence of such a path becomes more probable the higher the 
dimensionality of the moduli space is, but the danger of getting trapped in a meta-stable minimum 



e ~ ee where e = AVi/Vi; thus e 3> e is identical with our assumption of flat potentials. 
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remains. 

Within staggered inflation, we compute observables, such as n s and r, that are dominated by e; 
a comparison with the WMAP5 data [46 1 reveals that hilltop potentials are already ruled out at the 
2cr-level, while linear potentials are in excellent agreement with observations. We also re-examine 
briefly inflation driven by tachyons, as proposed in 211 ], and conclude that a constant decay rate 



is already ruled out in this setup while a serial condensation is still viable. 

We further comment on expected unique signatures of staggered inflation, such as additional 
gravitational waves and non-Gaussianities, that should arise whenever a field decays or stabilizes. 
To compute these signals, one has to go beyond the analytic formalism of this paper and study the 
decay of fields in detail. This would also provide a numerical check of the analytic framework. 

We conclude that staggered inflation offers a novel, compelling possibility to implement inflation 
within string theory, enabling the use of potentials that would be considered unsuitable otherwise. 

The concrete outline of this paper is as follows: we start by reviewing the background evolution 
of staggered inflaton in Sec. HU which is used in our discussion of scalar and tensor perturbations 
in Sec. IIII Al and IIIIBi We follow with an examination of concrete models, linear potentials in 
Sec. IIV Al and hilltop potentials in Sec. IIV Bl allowing for different slopes and masses in Sec. IIV ATI 
and IIV B~T1 The predictions are compared with WMAP5 constraints in Sec. IIV C( where we also 
re-examine inflation driven by multiple tachyons. We elaborate on the theoretical motivation for 
these models in Sec. |V] and point out open questions within staggered inflation in Sec. IVIl before 
concluding in Sec. IVIIl 



II. BACKGROUND EVOLUTION 

Consider M ~ O(10 3 ) uncoupled inflaton fields ifi, i = 1 . . .J\f, with potentials 

Vi = V - fi{w) , (1) 

where fi((fii) is a monotonically decreasing function, and the total potential is W = YJj V% 4 . We 
focus on those fields with a large Vq contributing the most to inflation, and ignore fields with 
much smaller Vq. Therefore, we set Vq identical for the fields we consider. Regarding the free 
functions fi(tpi), we are primarily interested in two cases: first, we identify inflatons with fields 
on the landscape which, by chance, encounter a flat region in the potential; we can then expand 
around such a position and get linear potentials fi = Ciifi. If, on the other hand, the fields are 

4 W should not be confused with the super-potential in SUGRA. 
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a) Potential needed for single-field inflation 
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FIG. 1: Schematics of potentials, shaded regions allow for slow-roll inflation: a) Potential needed for single- 
field inflation (such as the KKLMMT proposal) that is flat for Aip ~ O(10). b) A more generic potential 
on the landscape, with valleys, hills, and steep as well as some shallow regions, c) A region that is shallow 
until fend < 1 at which point the potential drops. If a field encounters this cliff, it decays (stabilizes) and 
p r is produced. We expand around this flat stretch, Vi rs Vo — c%^>% ■ d) A maximum, around which we also 
expand, Vi « Vo — mfipj/2. 

close to a maximum, we expand again and get quadratic potentials fi = (see Fig. [J for a 

schematic). The latter case appears less likely given a generic position on the landscape, but since 
no canonical measure is known (see Sec|V]for more details), we discuss both cases. 

In either setup, only narrow distributions of Cj or rrii are of interest to us, because in the presence 
of wildly different slopes, fields roll down quickly along the steeper ones and become irrelevant for 
inflation. For this reason we take fi = fj = f first, but relax this simplifying assumption as we 
proceed. 

We expect an expansion around either a maximum or a shallow region to be valid only for a 
limited stretch in field space before the potential changes drastically. To model this we use ([I]) 
only up until t/?j = <p en d -C 1; once a field equals (fend, we assume that it encounters a sharp drop 
in the potential, falls down and converts its entire potential energy into some other form p r , such 
as radiation. Since we would like to model flat stretches we demand 

- _ fi^end) 



«1. 



(2) 



In the following, we denote equality up to first order in small parameters with 
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If the fields start out from different initial values they encounter the drop in Vi one after the 
other and decay in a staggered fashion. If there are many fields so that several fields decay in 
a given Hubble time, we can apply the analytic formalism of staggered inflation as proposed in 



27, 



281 ] . We need to amend this formalism slightly to accommodate more general initial conditions 
as well as varying slopes and masses (Sec. IIV A 1 1 and IIV BT1) . 

Thus, let us briefly review (and extend) this formalism: first, we promote the number of fields 
M to a smooth, time dependent function, so that we can introduce the continuous decay rate 

r = -#>o. (3) 

This rate is not a free parameter but follows once the initial values for the inflatons are specified 
and a solution to their slow roll equations of motion is found. If inflation is to last several e- folds, 
we need 

£N = 1h ^ 1 ' (4) 

otherwise the majority of fields encounter (p en d in a single Hubble time. In the absence of a pre- 
inflationary model dictating the initial conditions, we distribute the fields evenly over the interval 
ip in i . . . (p end , namely 

/r\\ i Vend ~ fini , . lN / c s 

<Pi{p) = <Pini + — — (*-i); (5) 

this seems to be the most natural choice to us, but other initial distributions are possible (see 
SecJVjfor heuristic arguments motivating this choice). Note that quantum mechanical fluctuations 
generically introduce 5f® M ~ H. Hence, for hilltop potentials one should take ipi n { > H in order 
for the classical evolution to dominate 5 . 

Given these initial conditions, the fields evolve according to their Klein-Gordon equations. 
During slow roll, that is if the slow roll parameters e< = (dVi/dipi) 2 /(2W 2 ), e = £V well as 
Vi = {d 2 Vi/ dip 2 ) /W and r\ = ^ i rji are small, the equations of motion become 

3ff&~-|^. (6) 



Here, a dot denotes a derivative with respect to cosmic time and we set the reduced Planck mass 
equal to one, rn~ 2 = 8irG = 1. Eq.© has to be solved in conjunction with the Friedman equation 

3# 2 = Ptotal , (7) 



The COBE normalization of the power-spectrum implies ip^i > 10 6 in the present study, see section MI Al 
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where ptotai = Yli Pv> + Pr = Pi + Pr- The energy transfer from pi to p r manifests itself in the 

n n 

continuity equation, which becomes (see [27J and also [471 ] for related work where inflation is driven 
by a dynamically relaxing cosmological constant.) 



pi = -3H(pj+ Pl ) +MV (1 + 0(e)) . 



(8) 



Since V^T/*? , = has to hold, we necessarily need the additional component of the energy mo- 
mentum tensor to obey 



Pr 



-3H(p r + Pr )-MV (l + O(i)) 



In order for inflation to last several e- folds we need p r <C ptotai, that is 

«1, 



3/-, N Pr 

2 Pr+ PI 



0) 



(10) 



where w r is the equation of state parameter of the additional component, p r = w r p r . With these 
definitions, one can further simplify ([8]) and (j9|) to [27 ] 



pi ~ -2H(em + e)pi , 



Pr 



-2H ( - (1 + w r )p r - em Pi ) ^ 2H(em - e)pi 



(11) 
(12) 



where we ignored contributions that are second order in small parameters. As a result, p r ap- 
proaches a scaling solution that tracks the inflationary energy p r ~ Ej^2pi/(3 + 3w r ) for which 



em ^e. 



Further, the Hubble slow evolution parameter becomes 



27] 



H 

IP 



E + £ . 



Note that £ appears alongside £ and can be more important than the latter. As mentioned in 



this type of staggered inflation is reminiscent to warm inflation 
problems. 



(13) 



(14) 



351 ]. without sharing many of its 



III. PERTURBATIONS 



Scalar perturbations are discussed in [271 ] : here we summarize the main results of [27j, before 
proceeding to compute gravitational waves. Let's first briefly scrutinize the models of interest: we 
would like to focus on setups that possess a shallow potential, either because the fields are near a 



8 



maximum, or because such fields have, by chance, encountered a shallow stretch. We parameterized 
this requirement by demanding that e from is small, e.g. e ~ 10~ 2 . For the models of interest 
the slow roll parameter becomes of order e ~ ee <C s 6 . Hence, the contributions due to slow roll to 
observable parameters are negligible compared to those originating from the decaying fields, that 
is compared to those caused by £jv" ~ e if e is small. Thus, we drop all terms proportional to e. 



A. Scalar Perturbations 



As argued in [27[, isocurvature perturbations are suppressed in the present framework 7 , hence 
we focus on adiabatic perturbations; these are properly described by (the Fourier modes of) the 



Mukhanov variable v k = with z = 1/6, 6 2 = l/(3a 2 (l + w)) and w = Ptotal / ptotal |48l j. Here and 
in the following we use \c s \ — 1 (c s is the adiabatic sound speed), since p r oc pi during inflation. If 
i from (|14p evolves slowly, which is the case during inflation, we can approximate 

a(r) oc (-r) -(1+fi) , (15) 

where r = — oo ... is conformal time, adr = dt and d( )/dr = ( )', solve the equation of motion 
for v k 

< +(k 2 --)v k = (16) 
analytically in terms of Hankel functions 



48 



491 ] , match to the Bunch Davis vacuum in the far past 
(v k = exp(— ikr) /y/2k), expand the solution on large scales and compute the curvature perturbation 
on uniform density surfaces Ck- After some algebra, the resulting power spectrum for £ becomes, 
to leading order in e (see for details) 

n - ^irf (it) 



1 H 2 

i-rr 2 e m 2 z 



(18) 



where we restored the reduced Planck mass 8 . Note that the amplitude of the power spectrum is 
set by the COBE normalization P ? = (2.41 ± 0.11) x 1(T 9 [46]. With e ~ 10 -2 , this implies an 



6 For instance, in the linear case e ~ c 2 /(2A/"V 2 ), while e ~ ~ c/(2tp en dAfVo) from (|35p and e = c<p enc i/Vo, so 
that e ~ ee. Similarity, using (|53[) one can show that e < ee/2 for hilltop potentials. 

7 There is a possibility that isocurvature perturbations are produced during the short intervals when fields decay; 
we assume that this does not occur, but this caveat clearly warrants further study. 

8 For comparison, the power spectrum in [23], where fields start out from identical initial values and the slow roll 
contribution is kept, reads Vc, — H /(8tt m pl (£'y +e)), where 7 is of order one and it is set by F and the potential 
(see [13] for details). 
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inflationary scale of H ~ 10 6 , which in turn gives (fi n i ~ 5ipf M ~ H ~ 10 6 in hilltop potentials 



Further, the scalar spectral index becomes 



27] 



dlnP c 



"■ ~ 1 " d¥F (19) 

- (i-3)e, (20) 

where 

(21) 

TH 

measures the time dependence of the decay rate T, which also needs to be evaluated at around 
N ~ 60 e-folds before the end of inflation. Note that the condition for the validity of our approach 
(several fields ought to drop out in a given Hubble time for the smoothing of M(t) to be a good 
approximation) needs to be satisfied only when the Fourier modes C, k of interest leave the Hubble 
radius around N ~ 60; the reason being that once these modes leave the horizon they freeze out 

B. Tensor Perturbations 

Gravitational waves are less influenced by staggered inflation simply because they are not di- 
rectly coupled to scalars. To see this, consider the amplitude h of the two possible polarization 



states of a 
Laplacian 



general tensor perturbation hu = /i(i)e-^' X ' ) (x), where e« are eigen-modes of the spatial 



49l | . Fourier modes of this amplitude obey 



k 2 

h k + 3Hh k + —h k = 0; (22) 
a z 

hence, h does not couple to the scalar sector, but is determined entirely by background quantities, 
once initial conditions are specified. We can eliminate the friction term by switching to conformal 
time and introducing u k = ah k /2 so that 

ul+(k 2 --)u k = 0. (23) 



a 

Imposing quantum mechanical initial conditions in the far past u k = ex.p(—ikT)/y/2k, just as for 
scalar perturbations, and using fjl5f) . we find 



u k = ^Hll\-kr), (24) 
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where H$ is the Hankel function of the first kind with index v T = 3/2 + e. On large scales, this 
can be expanded to 

l u *l = \r- f W -t . ( 25 ) 



53] 



2vr \ 2 r 

where T is the Gamma-function. Thus, the amplitude of the tensor power spectrum becomes 



9 2~ 



(27) 



where we restored the reduced Planck mass in the last step and used T(3/2) = \/tt/2 9 . Further, 
the tensor spectral index reads 

dlnV T 

~ -2e. (29) 

Computing the ratio of ([27|) to (fT8|) yields the tensor to scalar ratio 

r - ^ (30) 
~ 16e. (31) 

Note that we still have r ~ — 8n^, just as for plain slow roll inflation, but both parameters are now 
determined by e — £j\f = T/{2H) instead of e. 



IV. MODELS 



In string theory, many degrees of freedom (branes, fluxes, size and shape moduli etc.) are present 
in the early universe, all of which can be dynamical. However, in most approaches to inflation all 
but one or two degrees of freedom are already stabilized during inflation. The argument for this 
approach is simple: even if all degrees of freedom are dynamical at the onset of inflation, only those 
acted upon by very weak forces, that is those with a very flat potential, will remain dynamical for 
a long period of time, hopefully driving inflation, whereas all the other ones quickly relax to stable 
configurations. Thus, by focusing on only one or two degrees of freedom, one assumes that flat 
directions are so rare that one would be lucky to encounter one suitable field for inflation. 

9 The additional factor of 2 in the power spectrum accounts for the two polarization modes of gravity waves. 



11 



However, this line of reasoning misses one crucial point: in the presence of many fields, each field 
sees the driving force of its own potential, but it is slowed down by the combined Hubble friction 
of all fields (this is the idea behind assisted inflation Qj). It is then not far fetched to employ many 
fields which, by chance, find themselves on a flat stretch on the landscape and can drive inflation. 
Of course, once they encounter a sharp edge, which should happen for Acpi <C 1, they rapidly fall 
down, convert some energy to other degrees of freedom such as radiation, and become irrelevant 
for inflation. For this type of inflation, one might expand around the flat stretch and use a linear 
approximation to the potential up until the cliff at •~p en d < 1, as indicated in section [TIJ In the 
next subsection, we discuss this model, followed by hilltop inflation. 



A. Linear Potentials: Vj, =Vq — api 

Consider N fields, with potentials Vi = Vq — ofi, that are valid up until (p en d, where we assume 
that the potential drops to zero. We distribute these fields according to (|5|) with (pint = 0, and 
apply the formalism of the previous sections to compute observables. To do so, we first need to 
extract T (we determine T without smoothing out the number of fields). Given the linear potential 
and a constant N ~ 10 3 we can compute the evolution of the fields from © to 

^^ >+ 7mt< (32) 

(whenever a field decays, N decreases by one in this expression). Thus, the time it takes for the 
field nearest to the edge to encounter the drop is given by ipj\f(At) = fend and reads to leading 
order in 1/N 



At ^^V^r (33) 



Consequently, the initial decay rate is 



r = ~ J- ~ - ( 34) 

N NAt </wv / 3W 



so that 



ftw = s^wW (35) 

Similarly, 5 can be computed using t = T 2 -N/N with N/N ~ l/(2AA(At) 2 ) ~ T 2 /2, so that 

YH f 1 

TH 2H 2 e 2 
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Further, the number of e-folds N ~ 60 becomes 

"end fend N 



N = I Hdt~- I 52^dw (37) 

J ini J ini • i "i 

~ , (39) 

where we used the slow roll approximation in the first line, the initial conditions from £< 1 
in the second one, and the large M limit in the last one. We can rewrite this expression as 
N ~ ip en dN / (2e) . Thus, with about a thousand fields M ~ 10 3 and flat potentials e ~ 10~ 2 we 
need the sharp drop to be around ip en d — 2eN/Af ~3.5xlCT 2 ^;lto yield around sixty e-folds of 
inflation. This is consistent with our expansion, since potentials should be well behaved over such 
small stretches in field space. Further, our requirement of a drastic change in the potential (such 
as a sharp drop) once fields get to (p en d appears natural. Given ([39]) . we can write 

e ~ — . (40) 
4N V ; 



The resulting scalar spectral index (|20|) . tensor spectral index (|29|) and tensor to scalar ratio (|3T[) 
are 

n.-l*~, (41) 

(42) 

r ~ 1 . (43) 

Note that r = — 4(n s — 1), just like in a single-field model with V = m 2 <f 2 /2, but the values for r 
and n s — 1 are smaller. A plot of the r-n s plane including the WMAP5 constraints [46J is in Fig. [2 



1. Different Slopes 

In realistic scenarios, it is unlikely that all potentials are identical, but we expect some spread 
of slopes. However, since fields with much larger slopes drop down quickly, we consider a narrow 
spread of the q's only; specifically, we take 

Ci = c(l + ljj) , (44) 

with I <C 1. If we use the same initial conditions as in ([5]) with = and use the above spread 
of Cj's, we implicity assume that fields with a larger slope are already further down their potential. 
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This seems reasonable to us, but the exact pairing of initial conditions and slopes could be chosen 
differently. Based on our choice, the number of e-folds becomes 

where we expanded for small /. Further, e — (1 + 2Z/3)/(4iV) and 6 ~ — 1, so that 

-- 1=! 4( 1+ t)' (46) 

B ^-^(l + §), (47) 
4 / 2A 

Thus, by increasing I, that is by spreading out the Cj's, the scalar spectral index and tensor to 
scalar ratio shift up along the r = — 4(n s — 1) line in Fig. [21 closer to the canonical chaotic inflation 
case. 



B. Hilltop Inflation: Vi = V Q - m 2 ip 2 i /2 

Consider a large number of fields in the vicinity of the maximum in their respective potentials. 
If we expand around a maximum, we have to be careful to avoid initial values too close to the 
top: fluctuations dominate over the classical evolution of the fields if fields are too close to the 
maximum. Since 5ipf M ~ H ~ we take (fi n i = <p en d/ot with a = 3.5 x 10 -4 (anticipating 

fend ~ 3.5 x 1CT 2 from below) in ©, that is 

W (0) = ^(l + (a-l)i^) . (49) 



As a result, cp%(0) - y»-i(0) ~ (Pend/N ~ 3.5 x 1CT 5 > </3j n j where M ~ 10^, as it should. Given the 
potential and initial conditions, the slow roll equations of motion integrate (for M ~ const) to 



2 



<ft(t) ~ ^(0)ev / ^ . (50) 



Using this result, we can, once again, estimate the decay rate to 



,2 



r = -^~^ Tl 2L_ (si) 



where 



t_ m2 in \ v l + ( a _l)(l_l/AA) y ' ~ a m 2 \ M ^ } 



14 



0.4: 



Staggered Inflation 



0.2 



0.1 



0.0 



0.4: 



0.3 : 



0.2 



0.1 



N= 50 60 



chaotic m 2 (|) 2 o 
staggered linear • 
staggered_hilltop o 

HZ 




Staggered Inflation from Tachyons 



N= 50 60 



chaotic m 2 4> 2 o 
constant decay rate • 
serial condensation o 




FIG. 2: We compare the predictions of staggered inflation with the WMAP5 constraints (reproduced from 
by permission of the AAS) and the canonical m 2 (p 2 chaotic inflation model. Shaded areas correspond to 
la and 2a confidence levels. Top panel: models explored in this paper with linear potentials (Vi = Vq — ctpi, 
section HV Ap or hilltop potentials (Vi = Vo — m 2 (p 2 /2, section lTVBP ; allowing for a spread of slopes or masses 
(I 7^ 0), shifts predictions up along the dotted lines. The hilltop case is borderline ruled out at the 2a level. 
Bottom panel: staggered inflation from multiple tachyons near a hilltop that condense due to fluctuations 



(not due to slow roll) as proposed in 2lJ, |27| , either with a constant decay rate (r = constant) or in a serial 
fashion (T = (t — l)" 1 so that s = 1/(3N)). The constant decay rate is already ruled out at the 2a level, 
while serial condenstaion is still viable. 

is the time it takes for the AA'th field to encounter the drop; here, we only kept the leading order 
in 1/Af in the last step. Consequently, 

m 2 a 

£ ^2W—1- (53) 
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Similarly, after some algebra, 5 = TH/ (TH) becomes 

S w -3 + - , (54) 
a 

to leading order in 1/M. In what follows, we would like to express e in terms of the number of 
e-folds N. The latter becomes 

pend N u y \f 

Jini i=l v i 

where m 2 /Vq = 2e/ip 2 nd and 

* = ln(«A0 " ln(a - 1) - ^ In rW(a _ 1)} J ■ ( 56 ) 

Generically, g is of order one for the models we are interested in and does not depend much on JV, 
for instance we get g ~ 1.0059 for J\f ~ 10 3 and a = 3.5 x 10 4 . As a consequence, we need again 
(fend ~ 2Ne/M ~ 3.5 x 10~ 2 to achieve sixty e-folds of inflation. Thus, we have 

all 

e ~ o « (57) 

y a-l2N 2N y J 

The resulting scalar spectral index (|20h . tensor spectral index (|29j) and tensor to scalar ratio (|3ip 
are 

" = ^-5' (59) 



Plotting the above in the r-n s plane including the WMAP5 constraints 
this model is already borderline ruled out at the 2cr-level. 



46J in Fig. [2] reveals that 



1. Different Masses 

Similar to Sec. lIV ATI we would like to allow for slightly different potentials for each field, that is, 
slightly different masses. Choosing m 2 = m 2 (l + li/N) with I <C 1 in conjunction with unchanged 
initial conditions, we implicitly assume that heavier fields are further down in the potential, which 
seems reasonable to us. Just as in the linear case, one can show that 5 ~ —3 remains unchanged to 
leading order in l/M. Thus, predictions stay on the same r(n s ) line in Fig. [2J Further, e becomes 

" '" 2 (l + O-^r, (61) 



2V N K a -I 
16 



and the number of e- folds reads 



N = ^fg(a,AT,l), (62) 



where 



~ 9 =h^ x T^ ln {w))- (63) 

The latter expression can be expanded to g « 1.0059-0.2510/ p» 1-1/4 for N = 1000, a = 3.5 xlO 4 
and Z < 1. Thus e « (2JV) _1 (1 + 3Z/4) and as a consequence, predictions move further away from 
the desired region in Fig. [2j 

n s -l«4(l + 4), (64) 
-T--^fl + /f)> (65) 



C. Discussion 

The predictions made in (|4ip - (|43p and (|58p -(|60p are generic for staggered inflation: the main 
assumption was that a large number of fields, J\f ~ 0(1O 3 ), encounters a flat, sub-Planckian 
stretch in field space before the potential drops off again. Flat means that the standard slow 
roll parameters of assisted inflation and e = AVi/Vq are small 10 . We expanded the potentials, 
keeping only the linear or quadratic terms, depending on whether the fields are close to maxima 
or generic flat stretches. Since only fields with similar shallow potentials contribute to inflation 
for an extended period, we focused on fields whose potentials are similar. After distributing the 
fields evenly, which appears to be a natural starting point, we observed that the slow roll dynamics 
causes fields to decay during inflation, one after the other; this determines e as well as 5, which are 
not free parameters. The scalar and tensor spectral indices, as well as the tensor to scalar ratio 
follow directly and can be expressed in terms of the number of e-folds N ~ 60. 

The hilltop case is already borderline ruled out at the 2a level according to the WMAP5 



constraints 



461 ]. see Fig. [2] top panel. If we allow for a slight spread of masses as in Sec. IIVB 11 



while retaining the same initial conditions, the predicted values of r and n s move up along the 



10 The latter parameter guarantees that e <C e, so that effects due to staggered inflation dominate; otherwise, one 
would have a standard assisted inflation setup. 
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r(n s )-line, further away from the desired region. Thus, the only way to reconcile a hilltop model 
with observations in staggered inflation relies in fine tuned initial conditions, quite different from 
the even distribution we used. The more generic linear case (see Sec. |V} is still viable. For 
comparison, we plot a single-field model with m 2 ip 2 potential. Note that both, staggered linear 
inflation and chaotic m 2 ip 2 inflation, have r = — 4(n s — 1), but the magnitude of r and \n s — 1| 
is smaller in staggered inflation. However, if we allow for different slopes as in Sec. IIV A 1\ the 
predictions move up along the r(n s )-line, closer to the canonical chaotic inflation case. Thus, the 
two models become degenerate in this plot, but could be discriminated by an observation of nx- 
An advantage of staggered inflation with the chosen parameters over such a large-field model is 

0, Q,Q, 



the possible avoidance of the ry-problem 



since no field has to traverse a distance longer 



than (fend ~ 3.5 x 10 2 ^ 1. If we either increase the number of fields, or we decrease e, <p en d 
becomes smaller, without affecting the above observables. 

1. Concrete Model: Inflation from Tachyons 

In the bottom panel of Fig. Owe plot n s and r for staggered inflation from multiple tachyons, a 
model proposed in j^ll. From a higher dimensional point of view, the tachyons are related to 
brane anti-brane pairs (Dp — Dp). By focussing on the Abelian part of the U(M) x U(N) gauge 



501 ] for a derivation as 



symmetry, Davis and Majumdar arrived at M tachyons with potential (see 
well as a general review on tachyons) 

W = Mt p - Cl N 2 + c 2 Yl l^l 4 + (H 6 ) • ( 67 ) 

1=1 8=1 

This potential is only valid in close proximity to (tpi, V>.\r) = 0. Here r p is a model dependent 
brane tension of order one while c\ ~ 0.87 and C2 ~ 0.21 [501 ] . The hundreds of tachyons start out 
close to the top of their potentials and they condense, that is decay, in a staggered fashion due to 
fluctuations, not because of slow roll. Thus, one can give all fields identical initial values without 
loss of generality. 

We arrive at the same expressions for n s , r and ny as in (|20|) . (|31|) and (|29p . if the slow roll 
contributions are negligible; disregarding these contributions is well motivated within this scenario, 
because the employed potential ceases 

we treat the decay rate as a free function, because we lack a good understanding of tachyon 



11 The tensor to scalar ratio reads r = 16(e7 2 + e) if slow roll contributions were included; here, 7 is a parameter of 
order one that can be determined from the potential as in [2^ ]. 



18 



condensation 12 ; two simple choices for T where proposed [2ll . |2jj: first, a constant decay rate 
r = const so that 5 = and e = 1/N; second, T = (t — t)" 1 where t = const leading to serial 
condensation so that 5 = —2 and e = 1/(3N). The constant decay rate is already ruled out 
at the 2a level, while serial condensation of tachyons provides a model in good agreement with 
observations Fig. [21 For further constraints on Af and r p , based on the requirement of achieving at 
least sixty efolds of inflation, we refer the interested reader to j^ . 



2. Summary 

To summarize, current observational bounds already rule out several models of staggered in- 
flation driven by fields close to maxima in their potentials, albeit a model with tachyons that 
condense serially is still in good agreement with observations. Also, the most generic case is still 
viable: many fields roll along flat stretches that can be approximated by linear potentials up until 
Aifi ~ 10~ 2 where a sharp drop follows. In the next section, we elaborate more on this model and 
speculate how it could arise in a realistic model on the landscape. 



V. STAGGERED INFLATION ON THE LANDSCAPE? 



The models studied so far appear to have ample freedom regarding the slope of potentials as 
well as initial conditions. If one has the freedom to choose either of them at will, one clearly looses 
any predictive power. In order to make the precise predictions of the last sections, we employ linear 
or hilltop potentials with a narrow distribution of slopes/masses as well as an even distribution of 
initial field values over the allowed interval. Here, we would like to re-examine how natural, or fine 
tuned, our preferences are. 

First of all, we need around a thousand fields for staggered inflation to work. Such a large 
number seems intimidating at first, but it is actually in the right order of magnitude for string 
theory, which contains many degrees of freedom (shape and size moduli, fluxes, branes, etc.). 
Instead of stabilizing all but one degree of freedom, as is usually done in attempts of implementing 
inflation into string theory , we let them evolve during inflation. This is similar in spirit 
to AA-flation [lj] (~ 1000 axions), inflation from multiple M5 branes (~ 100 fields) or chain 



39f| . Indeed, the mere fact that many degrees of freedom are involved in inflation appears 



inflation 

12 As a consequence, (p)reheating is presently not understood in this proposal and might indeed be problematic due 
to couplings to hidden sectors. 
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less fine tuned to us. 

How natural are the potentials we consider? We primarily assumed that there are some flat 
stretches over short intervals A(fi <C 1 followed by sharp drops with large slopes. Such potentials 
are less fine tuned than the ones employed in large-field models with a single inflaton field: in these 
setups potentials need to be flat over A(p ~ 10, which is hard to achieve, since quantum corrections 
usually lead to large contributions once Aip ~ 1 (the r] problem). Stabilizing all but one modulus 
and the necessity of keeping corrections under control are the reasons that single-field models like 
the KKLMMT proposal [5] appear rather baroque 13 . 

In contrast, we rely on large corrections to potentials that cause a given field to evolve fast once 
Aipi becomes of order 10 -2 . All we assume is that the subsequent evolution is rapid, leading to 
stabilization or decay of the field while releasing energy in the form of new degrees of freedom, 
e.g. radiation. Thus, we incorporate the notion of moduli trapping on the landscape 



52, 



but instead of assuming that the stabilization epoch ends before inflation commences, we propose 
that fields are still stabilizing during the last sixty e-folds. The completion of moduli stabilization 
goes hand in hand with the end of inflation and reheating, tying together two previously distinct 
phenomena. In this sense, staggered multi-field inflation is more economic than single-field models 
which require an additional (unobservable) earlier phase of moduli stabilization. 

Let's turn our attention to initial conditions; we employ an even distribution over regions where 
the potential is flat. How generic is this distribution? Addressing this question is challenging, due 
to our inability to assign a meaningful measure on the space of initial conditions (or the landscape 



m 



general for that matter; see i.e. 



54 



55, 



60H for a short review and 



5a ] for different proposals, 
61| for more comprehensive lecture notes). Nevertheless, we can give some plausible, heuristic 
arguments: consider a pre-inflationary phase where the potential energy is negligible compared to 
the kinetic energy of fields 14 . Such a regime arises naturally after a Hagedorn phase 15 , which has 



recently 



e.g. 



62, 



gained 



63, 



a lot of interest in String Gas Cosmology and bouncing models of the Universe (see 



641 ] for reviews). Since the fields move freely, they randomly scan all regions of the 



landscape accessible to them. While the universe is expanding (non inflationary), the kinetic energy 



13 These models have the advantage that all moduli are guaranteed to be securely fixed at late times, as they should 
be. See 0, 0] for estimates of corrections to the inflaton potential. 

14 If potentials are unbounded from above, they would introduce hard walls from which the fields get reflected, 
restricting the accessible regions on the landscape. There are models employing such features, such as the new- 
ekpyrotic scenario, see [d5| for a review. 

15 A phase where an increase of energy does not cause a corresponding increase in the temperature, but the excitations 
of long strings. This feature of a limiting Hagedorn temperature [661 ] is a generic feature of string theory, which 
encompasses T-duality. 
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redshifts until the potential energy becomes comparable to the kinetic one. At this instance, fields 
are at random positions on the landscape, a few already close to a minimum (or maximum), but 
the vast majority of fields find themselves on intermediate stretches that have non-zero slope. If the 
slope on such a stretch is large, the field relaxes quickly to a minimum (it gets stabilized /decays), 
but if the slope is flat enough, the evolution becomes over-damped and the field enters a slow roll 
phase. After a few Hubble times, the fields with the smallest slopes dominate the energy density of 
the universe and start driving assisted inflation [jj. Of course, only those with comparable (large) 
energy densities need to be considered here. Thus, we arrive dynamically at a setup similar to the 
one considered in this article, starting from a generic initial state. 

It seems plausible that intermediate stretches occupy more ground on the landscape than the 
very flat regions near extrema; hence, we expect to find the majority of fields in regions where linear 
potentials serve as good approximations. Thus, this line of reasoning favors the linear case over 
hilltop potentials 16 ; this is good news, since we saw that hilltop potentials are already ruled out at 
the 2<r-level. However, since we do not have a canonical measure on the landscape, these statements 
are mathematically ill defined. Therefore, we examined both cases in this paper. Nevertheless, if 
the fields have a periodic field domain, which is the case if the inflatons were for instance identified 
with axions as in AA-flation, one could make the above statements more concrete, since measures 
can be defined properly (we leave this to a future study). We anticipate that one is naturally led to 
an even distribution over the entire (finite) domain just before the fields start to see their potential. 
In this light, our choice of an even initial field-distribution is well motivated. 



VI. OUTLOOK 



Here, we summarize the main open questions regarding staggered inflation. 

a. Validity of the Analytic Framework: our predictions rely on taking a continuum limit that 
leads to a smooth decay rate T. This approximation becomes exact in the limit of M — > oo while 
keeping W and T fixed, but it remains to be seen how low M and T can be, before effects caused by 
the stabilization/decay of individual fields become observable. The physical reason for additional 
effects is the violation of slow roll by the decaying fields, causing i.e. a ringing in the scalar 



Classically, fields extremely close to a maximum would be the ones which survive the longest, and hence one might 
be tempted to conclude that those fields are the ones most likely to drive the last sixty e-folds of inflation. However, 
quantum mechanical fluctuations displace these fields so that they start rolling classically within a single Hubble 
time. They will then encounter cliffs on comparable time scales as fields on linear stretches. Thus, it is unlikely 
that a large number of fields simply hangs around near a maximum until late times. 
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power-spectrum, similar to the ringing after a sharp step in the potential of a single inflaton field 



in [26j] a simplified two field model was investigated numerically and indeed a ringing was 
found, but it turned out to be damped more than in the single-field case. If many fields decay in 
a given Hubble time, those features will be washed out even further, but exactly how many fields 
are needed is unknown. A numerical treatment should be performed, but a better understanding 
of the fields' decay is needed beforehand. 

b. The Nature of the Fields' Decay and Reheating: so far, an efficient and fast energy trans- 
fer from the decaying/stabilizing inflatons to radiation or other types of matter with a constant 
equation of state is assumed. The exact nature of this process is model dependent and differs 
considerably from case to case. It would be instructive to have a concrete toy model, enabling the 
numerical check of the analytic framework; the extensive literature on different reheating mech- 
anisms after inflation (see e.g. 



49J for a review) should be a valuable resource in this endeavor. 
Further, given a concrete model, the computation of additional contributions to obervables, i.e. 
non-Gaussianities and gravitational waves, becomes possible. 

But is reheating of standard model particles for the model at hand feasable in the first place? 
Due to the presence of a large number of fields, many of which are conceivably coupled to hidden 
sectors, one might fear that reheaing hidden sectors is a severe problem, just like in AA-flation 231 ] . 
Fortunately, there is a promising way of incorporating reheating into the framework of staggered 
inflation: since the fields decay one after the other, only the field(s) that decay last need to decay 
into (extensions of the) standard model degrees of freedom. Since we considered models on the 
landscape, we shoulc 
that were used in 



67, 



include all possible degrees of freedom, among which there will be i.e. those 



63, Isalzdl 

to drive inflation with the field content of the MSSM (see also 



711 ] for a multi- field model). If the last fie" 



freedom would be predominantly reheated 



d to decay is of this type, standad model 



69, 



degr ees of 



72J. The 



691 ] . for instance via instant preheating 
decay products of the previous fields (which are most likely in hidden sectors) would have been 
diluted to some degree at the end of inflation and might very well account for dark matter 18 . 



Another option would be to reheat via MSSM flat directions (see i.e. [73J and follow ups). 



17 This line of reasoning is one of the main motivations for implementing Inflation within the MSSM; note that in 
order to realize the low-scale inflation of i.e. [67l . [68l . [till [toI | one generically needs a prior high-scale inflationary 
era of the type we consider. Hence, low- and high-scale inflation complement each other. 

18 Depending on the concrete implementation of reheating standard model degrees of freedom, observables such as 
the scalar spectral index might differ from the ones predicted here: for instance, if a staggered inflation is followed 
by an extended low-scale, single-field inflationary phase driven [dT|. |68[ 1(391 [toI] that lasts sufficiently long, the power 
spectrum could be entirely determined by the latter. However, even in this case an understanding of the preceding 
staggered inflationary phase is useful, since it determines the initial state of the later low-scale inflationary stage. 
A study of reheating after multi-field inflation is in progress [33 ] . 
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c. 



Isocurvature Perturbations and Non-Gaussianities: within the analytic framework, isocur- 



vature perturbations are suppressed [271 ). leading to the expectation that non-Gaussianities are 
suppressed too [2^], just as in models of multi-field slow roll inflation 0,0]. Deviations from this 
expectation are feasible during the short instances when fields decay, leading to the production of 
non-Gaussianities similar to the ones generated in single-field models with a sharp step in the po- 



tential 



58 



59l | . or during (p)re-heating 



82 



83 



84 



85 



86, 



871 ] . Since fields are decaying/stabilizing 



throughout inflation, sizeable non-Gaussianities might result, which could very well dominate over 
primordial ones. 



d. Additional Gravitational Waves: it is known that gravitational waves can be p roducec 



sub-Hubble scales during (p)reheating (for a selection of numerical studies see 



3, 



80, 



57|, 



7_4j, 



76, 



77], 



on 



78], 



8l|). Because fields decay throughout inflation, we expect an additional gravitational wave 



811 ] ; since the 



spectrum resulting from the superposition of spectra similar to the ones found in 
process operates within the horizon, the peaks of the spectra are always generated at the same 
fraction of the Hubble radius whenever a field decays. Once produced, inflation stretches them to 
super Hubble scales, just as it does with gravity waves seeded by quantum fluctuations. Further, 
since each decay resembles the previous ones, the amplitude at which gravity waves are produced 
should not change either. As a consequence, the resulting superposition of power-spectra should 
lead to a spectrum with a spectral index quite different from ~ — 2s. Whether the primordial 
powerspectrum or the one from decaying fields dominates, dependents on the detailed physics of 
the fields' decay mechanism. Once a concrete model within string theory is constructed and the 
dacays are better understood, the additional contribution to Vt is computable and could serve as 
a smoking gun for staggered multi-field inflation. 

e. Feasibility of Staggered Inflation on the Landscape: the arguments within this paper on the 
occurrence of staggered inflation on the landscape were heuristic, due to the fact that a canonical 
measure on the landscape is still missing. Only after such a measure is known, one can properly 
estimate whether or not staggered inflation is generic, and whether or not initial conditions are 
fine tuned. Nevertheless, the heuristic arguments put forward in this paper are promising. 

/. Generalizations: so far, staggered inflation has been examined for separable potentials, 
canonical kinetic terms and a flat field-space metric only. String theory on the other hand permits 



more general scenarios; i.e. it wou 



DBI inflation 



89 



90 



91 



92 



93 



d be interesting to study the effects of decaying fields during 



94|] or k-inflation [95j, |9J] 
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VII. CONCLUSIONS 



Staggered multi- field inflation incorporates the notion of decaying/stabilizing fields during in- 
flation, a feature that has been largely ignored in the literature. We extended the formalism to 
incorporate more general initial conditions and potentials, enabling the discussion of generic multi- 
field models as expected in the landscape picture. The potentials we considered are flat over short 
intervals, before a sharp drop is encountered. Such potentials are less fine tuned than the one in 
single-field models, which require flatness over long distances in field space. After spreading out 
the fields evenly over the allowed interval and using either linear or hilltop potentials with some 
spread of slopes or masses, we computed observable parameters such as the scalar spectral index, 
the tensor to scalar ratio or the tensor spectral index. These parameters are dominated by the 
contributions due to staggered inflation e ~ T/(2H) oc Pradiation/ Ptotai > SslowrolV, e is not a free 
parameter but follows from the slow roll evolution. 

We found that hilltop potentials are already ruled out at the 2a level by the WMAP5 data, 
whereas linear potentials, which are better motivated, are in excellent agreement with observations. 
We also reconsidered staggered inflation driven by multiple tachyons [21], [2^], where T(t) is an 
unknown function. We find that constant decay rates are already ruled out at the 2a level, but a 
serial condensation of tachyons is still viable. 

One drawback of multi-field inflation is the apparent flexibility to choose initial values and 
potentials at will, increasing the tunability compared to single-field scenarios. As a result, multi- 
field models are in danger of loosing their predictability. To minimize fine tuning, we investigated 
which potentials and initial conditions appear to be generic on the landscape and gave arguments in 
favor of our choice. However, the absence of a canonical measure on the space of initial conditions 
and the landscape renders these arguments heuristic and open to debate. This uncertainty can be 
alleviated in certain models, for instance if the inflatons were axions with a finite domain. 

We concluded with a summary of open questions, ranging from a numerical validation of the 
analytic framework, to potentially new sources of gravitational waves and non-Gaussianities; the 
latter ones may very well serve as a smoking gun for staggered inflation, but require a better 
understanding of the fields' decay /stabilization mechanism before they can be computed. In this 
light, a more concrete implementation within string theory of this new type of inflation on the 
landscape is of great interest. 
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